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Abstract 



Bose-Einstein condensation (BEC) and evaporation of transverse photons 
from the Bose condensate is studied in the case when the density of plasma 
does not change. The generation of the longitudinal photons (photonikos) by 
the transverse photons (photons) in terms of the Cherenkov type of radiation 
in a uniform plasma is demonstrated. The Bogoliubov energy spectrum is de- 
rived for photonikos. A new physical phenomena of the " Compton" scattering 
type in nonlinear photon gas is discussed. To this end, a new version of the 
Pauli equation in the wavevector representation is derived. The formation of 
Bose-Einstein condensate and evaporation of photons from the condensate is 
investigated by means of the newly derived Fokker-Planck equation for pho- 
tonikos. The relevance of this work to recent discovery of black hole X-ray 
jets is pointed out. 



Submitted to Physical Review A 



1 



Typeset using REVTgX 



I. INTRODUCTION 



An immense amount of research has been carried out on propagation of relativistically 
intense electromagnetic (EM) waves into an isotropic plasma demonstrating that the rela- 
tivistic oscillatory motion of electrons causes a whole set of interesting and salient phenomena 
[1]- [6], relevant to the study of laser accelerators of electrons, ions and photons, laser fusion, 
nonlinear optics, etc. Some of them have already been confirmed by experiments due to the 
recent progress in compact, high-power, short pulse laser technology. In our previous paper 
[7] we have suggested a new mechanism for ultrahigh gradient electron and ion accelera- 
tion. In addition we have shown over lO^G generated magnetic fields by the nonpotential 
ponderomotivc force. 

The above treatments were restricted to the case of monochromatic EM waves. However, 
the interaction of relativistically intense radiation with a plasma leads to several kinds of 
instabilities and the initially coherent spectrum may eventually broaden, and naturally for 
ultrashort pulses the initial bandwidth is increasingly broad. Moreover, in astronomical 
plasmas there are a variety of sources of radiation, and in this case we speak about the 
average density of radiation from all the sources and their spectral distribution. Hence, 
the natural state of the strong radiation of EM field is with a broad spectrum. In order 
to study the interaction of spectrally broad and relativistically intense EM waves with a 
plasma, it was necessary to derive a general equation for the EM spectral intensity. Such 
an investigation was reported recently [8], [9], where the authors derived a general kinetic 
equation for the photons in a plasma. 

It should be emphasized that the radiation can be in two distinct states. Namely, one 
is when the total number of photons is not conserved. A good example is a black-body 
radiation. Another situation is when the total number of photons is conserved. Both these 
states have been studied in several aspects, mostly for the weak radiation. 

It is well known that the nature of EM waves in a vacuum is quit different from the 
one in a medium. In the vacuum EM wave exists only in motion, however the light can 
be stopped (wavevector k = 0, u = u{0) ^ 0) in different mediums and wave-guides. As 
was shown in Ref. [10] the photons acquire the rest mass and become one of the Bosons in 
plasmas and posses all characteristics of nonzero rest mass, i.e., we may say that the photon 
is the elementary particle of the optical field. 

Besides, a several reviews and books have been published on theory of Bose-Einstein con- 
densation (BEC) in a quantum Bose liquid [11] and in trapped gases [12], [13]. Kompaneets 
[14] has shown that the establishment of equilibrium between the photons and the electrons 
is possible through the Compton effect. In his consideration, since the free electron does not 
absorb and emit, but only scatters the photon, the total number of photons is conserved. 
Using the kinetic equation of Kompaneets, Zel'dovich and Levich [15] have shown that in 
the absence of absorption the photons undergo BEC. Such a possibility of the BEC occurs 
in the case, when the processes of change of energy and momentum in scattering dominates 
over the processes involving change of the photon number in their emission and absorption. 

Very recently it was shown that exists an another new mechanism of the creation of 
equilibrium state and Bose-Einstein condensation in a nonideal dense photon gas [16]. More 
importantly a new effect was predicted in the same paper, namely that the inhomogeneous 
dense photon gas can be found in the intermediate state. 
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In the present paper, we consider the BEC and evaporation of the transverse photons 
(photons) from the Bose condensate. In our study we assume that the intensity of radiation 
(strong and super-strong laser pulse, non-thermal equilibrium cosmic field radiation, etc.) 
is sufficiently large, so that the photon-photon interaction can become more likely than the 
photon-electron interaction. We will show that for certain conditions the variation of the 
plasma density can be neglected in comparison with the variation of the photon density. In 
such case the elementary excitations represent the longitudinal photons (photonikos) , for 
which we will derive the well known Bogoliubov's energy spectrum. 

The paper is organized as follows. First in Sec.II a basic equations, describing the rela- 
tivistic photon-plasma interactions, is presented. The problem of stability of the photon flow 
is discussed in Sec.III. The derivation of the Bogoliubov energy spectrum for the photonikos 
is given in the same section. Then in Sec. IV, we derive the Pauli kinetic equation from the 
Wigner-Moyal equation. Section V is devoted to BEC. The Fokker-Planck equation for pho- 
tonikos, by which we discuss the possibility of the creation of Bose-Einstein condensate and 
evaporation of the photons from the condensate, is obtained in the same section. Finally, a 
brief summary and discussion of our results are given in the last section. 

II. BASIC EQUATIONS 

If the intensity of photons is sufficiently large, then the photon-photon interaction can 
become more likely than the photon-plasma particle interaction. Under these conditions, 
which we shall assume to be realized, we may consider the medium to consist of two weakly 
interacting subsystems: the photon gas and the plasma, which slowly exchange energy 
between each other. In other words, relaxation in a photon-plasma system is then a two-stage 
process: firstly, the statistical equilibrium is established in each subsystem independently, 
with some average energies = KbT^ and Ep = KbTp, where Kb is the Boltzmann 
constant, is the characteristic "temperature" associated with the average kinetic energy 
of the photon, and Tp is the plasma temperature. The photon effective "temperature" , in 
general, will differ from the plasma temperature. Slower processes of the equalization of the 
photon and the plasma temperatures will take place afterwards. 

In the following we will show that under some conditions the photon-photon interactions 
dominate the photon-particle interactions. When the photon-photon interaction takes place, 
the phases of the waves are, in general, random functions of time. We need therefore not 
be interested in the phases and can average over them. In such a situation the perturbation 
state of the photon gas can be described in terms of the occupation number N{k,uj,f,t) of 
photons, and one can study how these numbers change due to the processes of interaction of 
the photons with each other, or with plasma electrons. Note that N{k,LO,r,t) is the slowly 
varying function in space and time. 

Recently, a new version of kinetic equation for the occupation number N{k, ui, f, t) of 
photons, for modes propagating with the wavevector k and the frequency u, at the position 
r and time t, was derived in Refs. [8], [9], [17]. It should be emphasized that the previ- 
ous derivations of kinetic equation were based on envelope equations, were restricted to 
nonrelativistic plasmas and neglected the time variation of the photon (or the wavepacket) 
mean frequency. In contrast, our derivation avoids the envelope approximation, is valid for 
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a fully relativistic plasma, and includes space and time correlations. Since this equation has 
appeared previously in the literature [8], [9], [17], it is presented here without derivation 

-N{k, u, r, t) + -{k- V)N{k, u, r, t) - ulsin- (V,^ • V,^ - -— ) • p ^ J ' = , (1) 



where ujp = y^^~^) P = ''^oe is the electron rest mass, and noe are the non- 

equilibrium and equilibrium densities of the electrons, respectively, and 7 is the relativistic 
gamma factor of the electrons, which can be expressed as 



TTo- 1+/?/ fduj N{k,u;,r,t) 



where 6 — ^^"^ 2 and fi is the Planck constant divided by 27i 



moenoeC 

Equation (1) is the generahzation of the Wigner-Moyal equation for the classical elec- 
tromagnetic (EM) field. We specifically note that from this equation follows conservation of 
the total number of photons, but not the momentum and energy of photons, i.e.. 



N 



where coefficient 2 denotes two possible polarization of the photons. Hence, the chemical 
potential, /j,^, of the photon gas is not zero. 

— * 

If the local frequency cu and wavevector k are related by the dispersion equation, cu — 
u!{k), then the occupation number is represented as N{k,u!,f,t) — 27rN{k,f,t)5{u! — u!{k)) 
and Eq.(l) reduces to 

|iV(«, r- *) + c'kV^^ - .J«nl(V, ■ Vj) ■ = . (4) 

ut 00 z oj 

In the geometric optics approximation {sinx ~ x), Eq.(l) reduces to the one-particle 
Liouville-Vlasov equation with an additional term 

':N(k, f, t) + V)NCk,^, r, t) - ^ (V,p ■ Vj - ^ ' ^^^^^ = . (5) 



From this kinetic equation we can obtain a set of fluid equations. To this end, we 
introduce the definitions of mean values of the density, velocity and Q = 7^ — 1 of photons: 

The photon mean velocity is defined by 

^ 2 /■ d^k fdujkc^ ^ . 

and 
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hup f <fk f du;N{k,u;,f,t) 

V — n'^^P / 7^r\3 / n ■ y° 

moen-oeC J [zny J In uj 

Having these definitions, we shall construct the transport equations by the usual way. 
Namely, multiplying Eq.(5) by 1, ^ and integrating over the entire spectral range, in k 
and UJ, we obtain equations of continuity, motion and Q of the photon gas: 

din 

——^ + divn^u^ = , (9) 
at 

^ + divQu, = , (11) 
where we have introduced the following notations in dimensional units 

t^ = -.V2/^/^^^%M, (12) 



(27r)3 7 2% u 



„ 2 f d^k f duj /kc^ ^\2,,,r ^ ^ ^ /-.^x 
^■-=3/ (2^05/ 2^(--''') ^'ft-.'-.*)=%r,. (13) 

Obviously, is the pressure of the photon gas in unit mass. For G we have the following 
expression 

huip Up r d^k f du N{k,uj,f,t) fiuip uJp , , 



moeC^ n-oe J (27r)^ i 27r uj^ moeC^ < uj^ > noe 

From Eq.(ll) follows that the term on the right-hand side (RHS) is much less than the 
first term on the left-hand side (LFS), if G/^y^ <S 1. This inequality allows us to neglect the 
term on RHS in Eq.(ll). With this assumption from Eqs.(9) and (11) follows the frozen- in 
condition 

— = const . (15) 

It is important to emphasize that in Eqs.(l), (4), (5) and (10) there are two forces of 
distinct nature, which can change the occupation number of photons. One force appears 
due to the redistribution of electrons in space, Vrie, and time, The other force arises 
by the variation of the shape of wavepacket. In other words, this force originates from the 
alteration of the average kinetic energy of the electron oscillating in a rapidly varying field 
of the EM waves, and is proportional to V7 and -^K 
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III. STABILITY OF THE PHOTON FLOW 



In order to study the problem of stability of a photon flow, we consider the propagation 
of small perturbations in a homogeneous photon flux- plasma medium. First, we will derive a 
relation between the variation of photon and plasma densities, from which we will establish 
the condition that allows us to neglect the variation of the plasma density in comparison 
with the variation of the photon density. 

To show this, we suppose the temperature of electrons to be nonrelativistic, i.e., -C 
moeC^7. In this case, the equations describing a state of the electrons are as follows 

% + (^e • = -eE - moec^^l , (16) 



^ + ci^^;-^ne = . 17 
at moe7 

We linearize Eqs.(16) and (17) with respect to perturbations, which are represented as 
Pe = 5pe; '^e = '^Oe + ^ud 7 = 7o + ^7 = 7o + whcrc the suffix denotes the constant 
equilibrium value, and 5pei Sue, and SQ are small variations in the wave. Taking into account 
the frozen-in condition (15), 5Q can be expressed by the variation of the photon density as 
5Q — Qq— — '^^^^5n^. After linearization, we will seek plane wave solutions proportional 
to expi{qf— Qi). 

Here we consider the range of low frequencies for which the inequalities Q < qc < uOp^ 
are fulfilled. In this case the photon flow can no longer excite the Langmiur plasma waves, 
and the contribution of perturbation of the electron density is rather small in comparison 
to perturbation of the photon density as 

5ne ^ 7o - 1 c^Q^ /^gx 
noe 27^ a;2^ uq^ 

Therefore, in our consideration the density of the plasma particles remain constant. In this 
case Eqs. (1) and (4) contain the variation of the relativistic 7 factor alone. 

In the following we shall demonstrate a new phenomena in the photon gas, which origi- 
nates from the variation of shape of the wavepacket. Obviously, any waves, which can arise 
in the photon gas-plasma medium in this case, are the proper waves with the weak decre- 
ment of the photon gas. It means that these waves are associated with the photon gas, but 
not with the plasma. 

To this end, we linearize Eq.(4) with respect to the perturbation, which is represented 

as 

N{k,u;{k),f,t) ^ No{k,uj{k)) + 6N{k,u;{k),f,t) . (19) 

The result is 
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or after summation we obtain 

where 57 =45Q = 4/(0-^^^^. 

After integration over all wavevectors k, from Eq.(21) we get the dispersion relation due 
to the relativistic selfmodulation 

1 + 4/5 /,3,riVo-g + .72) _ ^o-(^-^72)| _ (22) 

270 -^ ^ uj{k + q/2) u;(k-q/2) ^' ^'^cuik) ^ ' ^ ^ 

which has, in general, complex roots. 

This equation has been derived and predicted by us first for the monochromatic wave in 
Ref. [18], and then for the broad spectrum in Ref. [8]. 

In Eq.(22) No{k) is the occupation number in the equilibrium state and is represented 

as 



No{k) = no,{2^al)-'/'exp{-^-^^) . (23) 



This is a spectral Gaussian distribution, with the average wavevector and the spectral 
width (jfe. 

We can rewrite Eq.(22) in another form, taking into account a pole in the integral 



n-qu^O , (24) 



where u = < c. 

Using the well known relation 



lim — ^ — = (p— — t7T6(x) , (25) 

where p denotes the prescription that at the singularity x — the principal value is to be 
taken, Eq.(22) is rewritten in the form 

r ^.Kfi + ,72) ^ N,;(k - m) _ , „ ^6) 
27§ J u{k)^ u{k + q/2) uj{k-q/2) ! ^ ' 

— # — * 

In the first and the third integrals wc now replace the wavevector k + q/2 by k, and in the 
second and the forth integrals k — q/2 ^ k. Assuming that k — ko + x, \ q |, | X |"C| ko |, 
after integration we obtain the dispersion relation 

{Q-qUgy-qW^^-a^q'^ 2^lc^u;{ko)^ dxz ^x.=^{n-^,) ' ^ ^ 



where Ug = qUg = qUg cosQ, N^ixz) = I dXx I dXyNoix), Ve = c^^^^^j^, 

a — 2u(ka) ' ^"^^ velocity Vs = ca/3^^^^^^;^|^ can be treated as the sound velocity of 
longitudinal photons, similar to phonons in the quantum liquid at almost zero temperature 
[11]. 

We now first neglect the small imaginary term in Eq.(27), and examine in great detail 
the following equation 

{n-qUgr = q'{V,' + a\'-Vi}. (28) 
This equation has a complex solutions, when the inequality 

yj > + aV , (29) 

for q < ^^^-^/Vl — V^, holds. The unstable solution can be written as Q = flf + iflff, where 

Qf — qUg cos © (30) 

and the growth rate for the unstable modes with q < ~ is 



Q/f = q^Vi - - a^q^ . (31) 

This instability for the monochromatic relativistic EM waves was predicted in Ref. [18], 
and for the broad spectrum &X'Vs = ot = ^ was disclosed in Ref. [9]. 

If f2/ >| I, then the solution (30) clearly describes the emission of the longitudinal 
photons (photonikos) by the bunch of the transverse photons (photons) inside a resonance 
cone (cos© = 0,1 /qUg), similar to the well-known Cherenkov emission of EM waves by a 
charged particles moving in uniform medium with a velocity larger than the phase velocities 
of emitted waves. Also, here is an analogy with the Landau criterion about the creation of 
elementary excitations. 

Prom Eq.(31) follows that there exists two range of values of the wavevector q. Namely, 

2ui, fT?2 



for all q less than %v V| — V/, the Cherenkov type mechanism leads to the excitation of 



photonikos (large scale waves). For the other, q more than V| — V^, the frequency 
Q, — 0,1 -\- iQff is real. 

We now consider the case, when > V| or u"^ < joalc^. If we write Eq.(28) in the 
moving frame (Q — qu ^ Q), and multiply the both sides by the Planck constant h, then 
we obtain the well known Bogoliubov energy spectrum 



where s = hfl, p = hq, and rrieff = fiuj/c^. 

Equation (32) has the same form as the energy of the elementary excitations in the 
quantum liquid at zero temperature, which was derived first by Bogoliubov. 



Por small momentum (p < meffV = meff\JV^ — V|) from Eq.(32) follows e = pV ~ Vgp. 
This expression show us that the coefficient V is the velocity of "sound" in the photon gas. 
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similar to phonons in the quantum medium. For the large momentum {p ^ nieffV), the 
photoniko energy (32) tends to p^ /2m,,ff. i.e., the kinetic energy of an individual photoniko. 

Thus, we may conclude that the elementary excitations in our case correspond to pho- 
tonikos. 

The decrement, ImQ, — Qff, of photoniko is derived from Eq.(27), and the result is 

/mfi = -^F(-i^ye.p{-^^^} . (33) 

wnere - c ^^^^^y f - ^ 32 U{fco) J UfeJ a^' 

Clearly the expression (33) is valid for the weakly damped oscillations, that is ReVL 3> 
ImQ. We note that the decrement has a maximum near the frequency of the Cherenkov 
resonance fl ~ qu cos Q. 



IV. PAULI EQUATION FOR THE PHOTON GAS 

In the previous section we have derived the condition (18) allowing the situation, in 

which the density of plasma almost does not change, to exist. In this case from Eq.(l) 
follows that photons with different frequencies and wavevectors scatter on the photon bunch 
(wavepacket), and the equilibrium state may be established. This is a new phenomena of 
the "Compton" scattering type, i.e., the statistical equilibrium between the photons and the 
photon bunch will establish itself as a result of the scattering processes. 

In order to better exhibit this, we shall derive the Pauli equation in the wavevector 
representation. Note that the Wigner-Moyal-Tsintsadze equation (1) is convenient for a 
study of non-linear processes in the photon gas in the case of so-called weak turbulence 
states. Such states of the photon gas-plasma means that the energy of the quasi-particle 
(photoniko) is small compared to the photon energies. 

Usually the weak-turbulence states arise in a medium, when a small perturbation grows 
with the growth rate smaller than the frequency, ImQ ^ ReQ. In this case we can use a 
series expansions in powers of amplitude of the perturbations. 

In the case of weak-turbulence, we can define the distribution function N{k,uj,r,t) in any 
approximation for the amplitude of perturbation by iterating the Wigner-Moyal equation. 
Prom Eq.(l), we now derive an equation in wavevector space in the limit of the spatial 
homogeneity, which may be identified with the Pauli equation. To this end, we take the 
Fourier transform of Eq.(l) to obtain 

d ~* / ^Tcc^ \ ~* 
—N{k,u,q,n,t) jN{k,uj,q,n,t) 

2 [ d^ql [ dVLI „s . ^ d ^N{k,u;,q-qf,n-nf,t) 

where N(k,uj,q,Q,t) — J dfJ dtN(k,Oi),f,t)exp(^—t(q ■ r — ilt)), and 5p(q,Q,t) — 

J dfJ dt5p(f,t)exp(^—i(q ■ f— Qt)^ are the slowly varying functions in time. 

Let us consider two cases. First, in the case of g = and Q = 0, after changing variables 
—qf qi and —Vtl flf, we rewrite Eq.(34) in the form 
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dt = -"^ / (2^0^ y 2^^^^-^' — ^K^) — ■ 

Here we have used the identity 

sin|(g-Vg + fi^) =zsinhi(g.V,- + fi|^) =zL. (36) 

For the second case, when q^O and O 7^ (in this case d In N/ dt <^ fl), Eq. (34) can be 
written as 

fn - 2 /" cPg/ /" c?^^ . ^, ^ A^(fe, a;, g - g/, Q - Q/) 

= -^^y (2^y — M-?/,-^^/)^ — — ■ (37) 

By iterating this equation in the case of weak-turbulence we can write an equation for the 
occupation number N in any approximation in the amphtude of perturbations. Expanding 
the integrand about q = q/,Q = Q/ and keeping the first term with N{k, cu, q — q/, fl — fl/) — 
{2TT)^d{qf - q)d{n/ - n)N{k,cu), we obtain 

N{k, u, q, O) = -Up o - ^ - ■ ^^^^ 

Substituting Eq.(38) into Eq.(35) we get 

dN{k,u;) , f <fq f dn f \ Sp(q,n) \% N(k,u;) 

dt -"^^J i27rrJ 2nu;in-q-ur ^ ' ^ ' 

Since we are considering the weak-turbulence, for the perturbation the dispersion relation 
Q = Jl(g) is valid, i.e., 

6p{q, n) = 27rS{n - n{q))6p{q) . (40) 

Noting the relation 

d 1 

smh{y—)f{x) = -{f{x + y)- f{x - y)} (41) 

and 6p{—q,—fl) = 6p*{q,fl) (where the asterisk * depicts the Hermitian conjugate), we 
obtain from Eq.(39) 

dN{k,u;)_u^ f d^q \Sp{q)\^ fN{k + q,LU + n) N{k,uj) 



dt 4 y (27r)3^(A; + g/2) ^ uj{k + q) uj{k) 

X 



r (i:'q \dp{q)\^ r I\ [k + q. uj + \l) l\[k,uj)^ 
J (2nY u(k + u;{k + q) u;{k) ^ 

I + I 1 (42) 

n-q-uf VL + q-un' ^ ' 

n relal 

Use of Eq.(25) in Eq.(42) yields the following relation 



where ui — (*^+^/^)^ and use was made of dispersion relation u — uik). 



10 



dN{k) ^ f d^q 



9t E/ J^s[W±{k + q,k)N{k + q)-W±{k,k + q)N{k)] . (43) 
Here we have introduced the scattering rate 

W. = l f'^j^l m±?-u,). (44) 

^uj{k + q/2)uj{k) 

Equation (43) describes the three wave interaction, as iUustrated in Fig.l. First (a) 
diagram exhibits the absorption of the photoniko by the photon, whereas the second (b) 
diagram shows the emission of the photoniko by the photon. In other words, the pho- 
ton passing through the photon bunch (wavepacket) absorbs and emits photonikos, with 

— * — * 

frequencies = =f(u; — cof) and wavevectors q = =f(/c — kf). 

Integral in Eq.(43) is the elastic collision integral and describes the photon scattering 
process on the variation of shape of the photon bunch. 

It should be emphasized that the kinetic equation (43) is a new version of the Pauli 
equation in the wavevector representation [19]. Note that the equation type of (43) has been 
obtained by Pauli for a quantum system, whereas Eq.(43), derived for the dense photon gas, 
is pure classical. This equation indicates that the equilibrium of the photon gas is triggered 
by the perturbation 6p = 6{n/no'j), in particular by the perturbation in the wavepacket. 



V. BOSE-EINSTEIN CONDENSATION IN THE PHOTON GAS 



As we have mentioned in the Introduction, Zel'dovich and Levich have shown that 
the usual Compton scattering leads to the EEC in the nonrelativistic photon gas-plasma 
medium. 

In this section, we will show that exists an another new mechanism of BEC in a non- 
ideal dense photon gas. As was mentioned above, if the intensity of radiation is sufficiently 
large, then the photon-photon interaction can become more likely than the photon-electron 
interaction. 

It is well known from plasma physics that the transverse photon can decay into the 
transverse photon and a plasmon (electron Langmuir and ion sound plasmons). For this 
process to take place it is necessary to have a fluctuation of the plasma density. 

In our consideration, as was shown in Sec.III, the density of plasma remains constant. 
So that only one mechanism (relativistic effect), which can cause BEC, is the decay of the 
transverse photon into the transverse and the longitudinal (photoniko) photons. That is, in 
the weak-turbulence limit this mechanism is realized by three wave interaction, as we have 
discussed in Sec.IV. Namely, the photon (a;, k) generates photoniko and scatters on latter, 
with the frequency cuf and the wavevector kf. In this case, the wave number and frequency 
matching conditions are satisfied, i.e., k = ki + q, and uj = These processes continue 

as a cascade {uol = uj// + fl/, k/ = kll + ql^ etc.) till the wavevector of the photon becomes 
zero and the frequency uj = ajpf-j^^"^. This cascade leads to the BEC of the photons and the 
creation of the photoniko gas. When the density of latter is sufficiently small, the photonikos 
may be regarded as non-interacting with each other, however they can exchange energies and 
momentum with a bunch of the photons by scattering mechanism, and finally such process 
can lead to the equilibrium state of the photoniko gas, with the Bose distribution 
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where Q is the Bogohubov spectrum (32). 

We note here that since the total number of photonikos does not conserve, the chemical 
potential of this gas is zero. 

The result of BEC is that the ground state is filled by photons with the total rest energy 
Eq = Nq^^, where Nq is the total number of photons in the ground state. 

We now express explicitly the energy of the ground state, Eq, through the total num- 
ber, A^o; and the volume of the photon gas in two limits. First, we consider ultrarel- 
ativistic case, when 7 = ^/T+~Q ^ Q}!'^ = {huJp/moC^f^^{No/Noef/^{l/V^/^), where 
cJp = {Ane"^ Nq,; / rriQ^.y/'^ , Nq^ is the total number of electrons. 

The chemical potential of the photon gas is 

Note that in quantum Bose liquid at T = 0, the chemical potential increases with increase 
of Nq, whereas in the present case /j, decreases as shows Eq.(46). 

We now find the pressure in the condensate, Pc, as Pc — —dEo/dV, and write the 
equation of state for given Nq, the result is 

P^yV3 ^ const . (47) 

Next in the nonrelativistic limit, we have 7 ~ 1 + {l/2)(hiJp/moc'^){NQ/NQe){^/V^^'^). 
Then for the equation of state and the chemical potential we get the following expressions 

^ ^^^^^ ^4g^ 

and 

The second term in Eq.(49) is always less than one. 

We note that when the conditions for the formation of BEC arise, then the photon 
gas in plasmas can be far from the equilibrium. The kinetics of this phenomena is a very 
important issue not only in the study of the contemporary problems of short pulse laser- 
matter interaction, but also for understanding of the processes, which take place on cosmic 
objects in the presence of strong radiation of EM field. 

The kinetic equation (43) is the simplest nonequilibrium model in which we can expect 
BEC. To show this, we introduce the number density of the photons in the condensate as 
riQ — limy_^oo Then we can write N{k, t) — A.'K^5{k)nQ{t), which means that the photons 
with a zero wavevector are in the condensate. 

We next discuss an another situation in which the photon with the wavevector k is 
collided with the photoniko with q. For such processes exists the probability of the condition 

— * 

k -\- q — Q tohe satisfied for the photon and photoniko. Thus the photons absorbing the 
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photoniko will pass to the ground state. Therefore, the occupation number N{k + q), which 
is under integral in Eq.(43), can be written as 

N{k + q) = 4Ti^S{k + q)ni{t) + N{k + q)^^_^ , (50) 

where ni{t) is the density of pairs, for which the equality k + q = holds. 

Substituting Eq.(50) into Eq.(43), for the Bose- Einstein condensate density no{t) we 
obtain 

g^ = "''"-"°"> +2/. ^fJ^^N(k,,t)W(Mk'-,l, (51) 
ot To Jki^o {2TTy J {2ny 

where r^' = f ^W{q) = ^uj^ J ' %^"'") ' ^^^^ 

— * — * 

replaced k + q ^ kt. 

We specifically note here that the first and the third terms in Eq.(51) are the sources of 
production of the Bose-Einstein condensate, whereas the second term describes the evapo- 
ration of photons from the ground state. 

As we can see from Eq.(51), the stationary solution may exist, and it is 

no = ni H ^ — n2 , (52) 

< r > 

where < r >-'^ ^ / with r(^/)-i = / ^W{y, kf-q). 

Equation (52) has an obvious physical meaning, which is the following: the departure 
of photons from the Bose-Einstein condensate should, under steady-state conditions, be 
completely balanced by their arrival from other modes. 

The problem of BEC and evaporation of the Bose-Einstein condensate we can investigate 
by Fokker-Planck equation, which we shall derive from the Pauli equation. Here we suppose 
that the wavevector q and the frequency Q of photonikos are small in comparison with the 
k and uj. This allows us to use the following expansion in the integrand 

^ ^ ^ (9 an d'^WN 

W{k + q, k)N{k + q)^ W{k)N{k) + Q-^{WN) |,=o +^^^ l.=o +- ■ (53) 

Substituting this expression into Eq.(43), we obtain the Fokker-Planck equation for pho- 
tons, which describes the slow change of the occupation number in the wavevector space 



where 



A= / -0^q^W{q,k) , (55) 



D^J = I k) . (56) 
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The quantity A is the dynamic friction coefficient, whereas D^j is the diffusion tensor in the 

wavevector space. 

dD 

Introducing the following definition = ■> ^'^'^ noting that the expression on 

RHS in Eq.(54) is divergent in wavevector space, Eq.(54) can be rewritten in the form 

dN d ^ 

where jj = —B^N — ^ff- is the photon flux density in wavevector space. The fact that the 
flux should be zero allowed us to express B^ and in terms of one another. 
The equilibrium distribution function we choose to be Gaussian 



k 



2 



N{k) - const • e . (58) 
Substituting expression (58) into the equation J= 0, we obtain 

B.al = ^A;,A, ■ (59) 
Finally the transport equation of the photon gas takes the form 

dt dk^ 2 y al ^ dkj^ ■ ^ ' 

In order to solve Eq.(60), we consider a simple model representing the diffusion tensor 
as Dy — DqSij, with Dq = const. Such situation is realized when u ■ cos© <^ qc^/cu. With 
this assumption Eq.(60) reduces to 

^ = alikN) + ^VlN, (61) 
dt dk ^ 

where a = 

To discover the physical meaning of terms on RHS of Eq.(61), we consider them sepa- 
rately. First we neglect the diffusion term in Eq.(61), and assuming k{0, 0, k), we get for the 
new function / = kN 

The general solution of this equation is an arbitrary function of ko = ke"'^ (where ko is 
the initial value of the wavevector), i.e., /(/co) = f{ke"'*). This function is constant on the 
curves ko — /ce"*, while at i — > 00 the wavevector goes to zero (/c — > 0). The meaning of this 
is that the occupation number would tend to peak toward the origin as 

N(k, t) = T = const ■ e"* . (63) 
k 

Thus we conclude that the friction effect leads to BEG of the photon gas. 
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It should be emphasized that the exponential decay of the wavevector k = /cqc""* is due to 
the linearity of the Fokker-Planck equation. Note that we have neglected all nonlinear terms 
in deriving the Pauli equation. The presence of nonlinearity could saturate the exponential 
growth of N{k,t). In the above approximation (63) the condensate formation time, tc, is 
defined by the relation a^c ~ 1- 

— * 

Second, we suppose that W{q,k) is the even function with respect to q. In this case 
A — 0, and we have 

^ = ^^IN . (64) 

dt 2 ^ ^ ' 

Assuming that initially all photons are in the ground state with A; = 0, and the occupation 
number of the photons is N — 47r^no(5(A;), then the solution of Eq.(64) reads 

N{k, t) = noe'^'/{27iDQtY/^ . (65) 

In the course of time, the number of photons with A; = decreases as t^^^'^. The number 
of photons in the surrounding wavevector space rises correspondingly, and initially peaked 
at the origin is to be flatten out. Let us determine the mean square wavevector from origin 
at time t. Prom expression (65) we have 

< >= Dot . (66) 

Thus, y/< k^ > increases as the square root of time. 

Eqs. (65) and (66) manifest that the evaporation of the photons from the condensate 
takes place. 

We now derive a relation between the diffusion time, tr), and the time of the condensation. 
Prom Eqs. (63) and (65) we obtain 

2al ^ ^ k' 

ana id 



Do 2Do ' 

where (7^ = as follows from the Wigner function N(k,r) — Noexp[—^ ~ Note 
that To is the initial cross size of the photon bunch. 
Pinally, we arrive at the desired relation 

J- - k-rl . (67) 

From here it is evident that for the condensation and evaporation of photons, it is necessary 
that the following inequality to ^ tc is satisfied. This is to be expected, as in realistic 
astrophysical objects with a strong radiation, as well as in laboratory experiments with a 
laser radiation, the following condition k'^r^ » 1 is vahd. 



VI. SUMMARY AND DISCUSSIONS 

We have investigated a class of problems involving the interaction of spectrally broad 
and relativistically intense EM radiation with a plasma in the case when the photon-photon 
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interaction dominates the photon-particle interactions. We have presented a new concept 
of the estabhshment of equihbrium between the photon and the wavepacket of EM field 
(the dense photon bunch). This is a fully relativistic effect due to the strong EM radiation. 
We have established the condition under which the variation of the plasma density can 
be neglected in comparison with the variation of the photon density. In such case the 
elementary excitations represent the photonikos, for which we have derived the well known 
Bogoliubov energy spectrum. We have studied the BEC and evaporation of the photons 
from the Rose condensate in the case when the density of plasma does not change. To this 
end, from the Wigner-Moyal-Tsintsadze equation [8], [9], [17] we have derived a new version 
of the Pauh kinetic equation for the photon gas. For the case, when the wavevector and 
frequency of the photoniko is small in comparison with the wavevector and frequency of 
the photon, we have derived the Fokker-Planck equation for photonikos. We have presented 
a simple model, which exhibits the possibility of the creation of Bose-Einstein condensate 
and evaporation of photons from the condensate. We think that these processes can be 
detectable in next generation experiments with appropriate instrumentation. In fact, a 
number of experiments have been carried out in which plasmas are irradiated by laser beams 
with intensities up to 3 ■ lO'^^W/cm^. At such intensities the photon density is of the order of 
~ lO^^cm"^. For the plasma densities up to rig ~ 10^^, we should expect that the Bose- 
Einstein condensate becomes observable. The theory developed in this paper should also be 
the case for astrophysical objects, such as a black hole. In this connection we speculate that 
the recently observed radiation from the black hole may be attributed to the evaporated 
photons from the Bose condensate. That is initially all photons fall into Bose-Einstein 
condensate, and then after a certain time, as discussed in this paper, some photons undergo 
evaporation from the condensate. These processes may also explain the observable variation 
in radiation intensity, from being undetectable to one of the brightest sources on the sky. 
Note that these sources can turn off for decades, and new ones are always being found. In 
addition Universe is filled with a gas of photonikos. This gas may play the decisive role in the 
expansion of the Universe. Moreover, it may also be useful for explaining certain processes 
in supernovae explosion. Finally, the present theory may also find a valuable application in 
the future space technology, as well as in nonlinear optics. 
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FIGURES 

FIG. 1. a) Absorption of the photoniko by the photon 
b) Emission of the photoniko by the photon 
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Fig. la Tsintsadze 



Fig. lb Tsintsadze 



